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, The article provides full-analytic gravitational wave (GW) forms for eccentric nonspinning com- 

pact binaries of arbitrary mass ratio in the time Fourier domain. The semi-analytical property 
5_( ' of recent descriptions, i.e. the demand of inverting the higher-order Kepler equation numerically 

but keeping all other computations analytic, is avoided for the first time. 
, The article is a completion of a previous one (Tessmer and Schafer, Phys. Rev. D 82, 

■ 124064 (2010)) to second post-Newtonian (2PN) order in the harmonic GW amplitude and 

OO ' conservative orbital dynamics. A fully analytical inversion formula of the Kepler equation in 

harmonic coordinates is provided, as well as the analytic time Fourier expansion of trigonometric 
functions of the eccentric anomaly in terms of sines and cosines of the mean anomaly. Tail terms 
are not considered. 
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1 Introduction 



In our previous article [I] - it will be called "paper 1" from now onwards - we provided the Fourier 
domain gravitational wave (GW) forms including the first post-Newtonian (IPN) corrections to 
QQ , the leading order quadrupolar contribution and IPN corrections to the Newtonian accurate orbital 

CO ' dynamics in terms of tensor spherical harmonics. In this article, we complete the amplitude to 

, 2PN order in harmonic coordinates and we include 2PN orbital dynamics as well. We regard 

nonspinning compact binaries with arbitrary mass ratio in quasi-elliptical orbits, where the orbital 
dynamics is well-known. Earlier works provided the solution to the evolution equations in terms 
of a Keplerian parameterization, perturbed by their post-Newtonian corrections. The result is 
the so-called quasi-Keplerian parameterization (QKP), to be found, e.g. in [U[31ll], which will 
be one basic for this article. 2PN corrections to the GW amplitude in the limit of non-spinning 
binaries are also dealt with in [S] where the wave forms are provided with the help of a quasi- 
' Keplerian parameterization, and recently, including spin, in [6] where the results are given in 

terms of coordinate velocities. We like to improve the use of those wave forms as we provide 
them as pure series of harmonics in the time Fourier domain, as this is one step closer to data 
analysis investigations. Let us now briefly summarise why 2PN corrections and the inclusion of 
eccentric orbits are necessary. Current data analysis investigation showed that, even for circular 
inspirals, the inclusion of 2PN orbital dynamics is necessary to guarantee a successful detection of 
the GW signal [H [71 [S] . Authors of a recent publication [5] investigated certain equal- mass black 
hole binaries in the millisecond GW frequency range in rotating clusters and stated parameter 
estimation errors for the initial eccentricity of ~ 10^^ for LISA. Therefore we claim that to 
perform an eccentric GW analysis for LISA is demanded. The GW energy flux in terms of its 
harmonic constitutive parts, was given to IPN order, in the extreme mass-ratio regime, for the 
first time - to our best knowledge - in |10) . First attempts, however, for a time Fourier domain 
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(TFD) of eccentric h^^ and h^, the two polarizations of the far-zone GW field, have been made 
in Reference where the periastron advance has been incorporated by hand and the stationary 
phase approximation has come into use, in [12] and |13j , where a spectral analysis for steady-state 
binaries in the simple case of Newtonian motion and amplitude has been performed. 

Recently in |14| . the authors furnished Newtonian accurate TFD wave forms as they incorpo- 
rated the sine and cosine functions of the eccentric anomaly as Fourier-Bessel series. Taking as 
starting point /i+ and themselves, having evaluated all appearing scalars (such as, for example, 
the scalar product of the basic vectors of the plane of the sky with the orbital velocity vector) 
being evaluated explicitly, would lead to rather complicated relations of rotation angles connect- 
ing the orientation of the binary's orbital plane with respect to the position of the observer (see, 
e.g. [15], and - including spin precession - [16] for higher orders of h-^- and hx) and this is not 
well-suited for a systematic TFD representation at higher orders. 

The problem of nonspinning compact binaries implies the conservation of the direction of the 
orbital angular momentum h, thus, it qualifies for a representation with the help of tensor spherical 
harmonics |17l [T51 [T^ . Reference [SU] provided IPN accurate spherical tensor components for 
nonspinning compact binaries in quasi-elliptic orbits (and quasi-hyperbolic ones as well) that we 
are on to verify in this article. Irrespective of the observer's orientation, the form the GW signal 
can be expressed with the help of these spherical tensor components and applied to give the 
polarizations h+ and [H] afterwards with the help of constant rotations. 

We provide an analytical expression of the eccentric anomaly u in terms of the mean anomaly 
M, and do likewise for the sines and cosines of arbitrary integer multiples of u at full 2PN order. 
By means of this, we are able to express the spherical tensor components in terms of infinite series 
in harmonics of the elapsed time. To which order the infinite summations have to be driven will 
be investigated. Similarly to Reference [55], tail terms occurring through 2PN order will not be 
discussed. 

1.1 Organisation of the paper 

The paper is organized in the following way. Section [5] collects the orbital elements of the quasi- 
Keplerian presentation of the 2PN-accurate orbital motion in harmonic coordinates. Section [3] 
summarizes the transverse-traceless (TT) projection of the far-zone GW field up to 2PN order 
corrections to the leading-order quadrupole field. The time derivatives of the source multipole 
moments are computed using expressions for the particle accelerations known from the literature, 
and the scalar components are projected out with the help of TT tensors in Section |4l The vital 
components of the Fourier decomposition, here at 2PN accuracy, of the GW field are summarized 
in Section [S] and the Fourier domain multipoles, incorporating purely 2PN conservative dynamics, 
are given structurally in its Subsection [531 The full explicit expressions can be found in Appendix 

m 
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1.2 Notation 



An overview of the quantities which are in use in this paper is given in form of a table of their 
associated abbreviations and scalings. 



u: eccentric anomaly 

M: mean anomaly 

h: orbital vectorial angular momentum, h := |h| 

xPN: x**^ post Newtonian 

r radial separation of the binary components 

Er "radial eccentricity" 

ct "time eccentricity" 

60 "phase eccentricity" 

Gr semimajor axis 

4) orbital phase 

V true anomaly 

$ total phase elapsed from a periastron to periastron time 

Af .mean motion, Af := where P is the time from periastron to periastron 

K. periastron advance parameter. K. := ^^J"^ 

7] symmetric mass ratio, rj :~ ''^^Jfi'^ , ~ mi + m2 as the sum of the masses 

Sm the difference in the masses, dm ■= wi — m2 

\E\ absolute value of the scaled orbital binding energy 

G Newtons gravitational constant 

c speed of light 



Multiple indices are used according to Thorne and Blanchet [T71 . In the final form of our 
results, energies are scaled with /xc^, separations with momenta with fic and times with 
such that all quantities are dimensionless. Lower index "0" labels the value of a quantity at time 
i = to- 



Schematic prescription of our calculation 



For convenience, we show a flow diagram of our computation. Preliminary results have been 
surrounded by boxes. 
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far-zone field: TD's (=time derivatives) of STF multipole moments and Jl 



to compute time derivatives: use quasi-Keplerian parameterisation 
to deal with scalars only: transform to spherical tensor components /'™ and 5'''^ 

4 



TD's of J and S as functions of eccentric anomaly u 



4 

decompose and sort appropriately 



TD's of / " and S . sum of even and odd functions of ' 



convert (symbolically) 
even terms: X)m>o cosmw, 
odd terms: ^^^q s,„ sin mit, 
4 



TD's of /'™ and 5*'"': sum of trigonometries of mu over positive m 



convert (symbolically) 
cosTOu: 'Yli'j><d"1i cosjA^, 
sinmu: X]j>o i sinj.M, 
{m e IN) 



TD's of I and S . sum of trigonometries of jA4 over positive j 



Fourier transform: straightforwardly performed. 
The single steps are being detailled below. 

2 The binary orbital motion: 2PN accurate Quasi-Keplerian pa- 
rameterization 

To compute the time derivatives of the multipole moments which we do in harmonic coordinates, 
we take the orbital parameterization at 2PN from [3], 



r — ar (1 ~ Cr cosw) , 
M = Af {t — to) ~ u~ et sin u -f- [Fy-u {v — u) + Fy sin v] 

2 TT 

— {(j) - M =^ V + (e'^A^) sin2u 4- {e'^g44>) sin 3?; , 



(1) 
(2) 

(3) 



where v = 2 arctan 



1/2 



tan ■ 



The explicit 2PN accurate expressions for the orbital ele- 



ments and functions of the generalized quasi-Keplerian parameterization, in harmonic coordinates, 
read 



1 



, 2(2|^l) , , ^ 4(2|S|) 

^ + e'^^ (-7 + 77) + ^^-^Ye 



16 



{2\E\h^ 



■(-4 + 777) 



(4) 



2 = 1- 2\E\h^ ' '2 



e^(^|24-4r; + 5(-3 + r7) i2\E\h') 



Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



5 



8 

M^{2\E\) 



jeO + 148 77 + 2 - i2\E\h^) (80 - 45 77 + 4 77^ 



32 



(4-777) , 



^/^Jl + .«(-15 + 7,) + .^(2|^l) 



2 r 



192 



8 



= (-5 + 277) 



128 



555 + 30 77 



2 6|ii;| , 4 -3|i?|2 ((977 - 22)6? +977 -21) 



We/ 
^IPN 



We/ 



et^ = 1 - 2\E\h^ + e 



2 , 2(2|i?l) / 



8 + 877- (2|£;|/i2)(-17 + 777) 



^4(W 

8 



4 [ 

12 + 72 77 + 20 772 - 24^(21^17^ (-5 + 2 77) 



(2|£;|/i2)(ii2_ 4777 + 16772)- 



16 



{2\E\h^) 



(-4 + 777) 



24 



= (-5 + 277) , 



E, 



mE\f 

{2\E\f ( ^l-2\E\h^ 



(-5+277) 



8 I 

/l2 4 



77 (-15 + 77) L 



2 r 



{2\E\K^ 



-(-5 + 277) 



15 



[2\E\h^y 



r(-7 + 27;) 



8 l(2|^|/72)' 

(2|£;|)^ r (i-2|£;|/i2 



32 



(2|i?|/72)^ 



■77 (-1 + 377) , 



e^^^l-2\E\h^ + e^^^'^^ ^ ' ^'^'^"^^ 



24+ (2|£:|;72)(_;L5^^^) 



^4(^1^ <! - 40 + 3477 + 1877^- (2|^|/i2)(160-3l 77 + 3 772) 



16 
1 



-(-416 + 91 77+ 157/) I , 



where we introduced the abbreviation 



We. Vl-e?- 



(5) 

(6) 
(7) 



(8) 
(9) 
(10) 

(11) 
(12) 
(13) 



(14) 



(15) 



There are 2PN accurate relations connecting the three eccentricities e^, Ct and e^. These relations 
read 



et 



er< 1 



1 



4c4 (2|£;|/i2) 



16 + 287? 



-(-30 + 1277) VW^M + (36 - 1977 + 6?72)(2|£;|/i2) 



(16a) 
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1 + e — - — rj 



1 



2 ■' ■ 32 (2|£;|/i2) 

+{-'n + iW){mh^) 



160 + 35777 - 1577^ 



(16b) 



Note that the set (|i?|, et, (/)o) is completely sufficient to describe the initial data of the system in 
its plane perpendicular to h, thus, we eliminate h for the benefit of et in each expression with 
the help of Equation ([5]). We could also have worked with the QKP in ADM coordinates, but 
this would have implied a transformation of the multipole moments. Reference |22| provided this 
transformation for non-spinning compact binaries, coming from a change in the positions x and 
the velocities v at 2PN. Reference [53] supplied the associated spin-dependent terms to 2PN (as 
we count the spins of order 0{e^)). Working in harmonic coordinates makes life more easy as 
harmonic coordinates are a kind of "natural environment" for the GW prescription and, there, 
those transformations disappear. In the next section, we list the relevant GW moments and their 
time derivatives (thanks to the work of Blanchet, Damour and Iyer ^51 [^ . Reviews of this 
subject can be found in [TH] and [H]). 



3 The relevant multipoles and their time derivatives in harmonic 
coordinates 

The far-zone regime GW forms, regarding only the instantaneous (non-tail) parts, can be expressed 
in terms of the mass- and current- type multipole moments jl7l I21| , 



G 



c^R' 
c 



(N)E 

1=2 

81 



1-2 



X,M,-2 it-R/c)Nj 



(i + iy. 



(17) 



Here, Ai is a niulti- index with Ai = ai...a/, indices with round brackets are symmetrized over, viz. 
for example ^(ij) = ^ {Aij + Aji), and the following definitions came to use: 



R 



= NR, 

— T 

= Na,...Na,_2 



(18) 
(19) 
(20) 

(21) 



R denotes the distance from the observer to the binary. The quantity N with multiple indices is a 
tensor product of components Ni of the normal vector N (representing the line of sight from the 
observer to the center of mass of the binary) from Na-^ to Nai_2, and the superscript (l) denotes 
the Z*^ time derivative. In the current case of 2PN accurate orbital dynamics, we have to evaluate 
the GW amplitude. Equation P?|) . consistently to 2PN order relative to the leading term, which 
kicks in at e^. The mass- type multipoles relevant for the above equation read [22| 



/i STFy 
1 



42 c2 



(29 - 87r])v^ - (30 - 4877) 



Gm\ 



1 



1 



c4 V504 



(253 - 183577 + 3545772)7;-^ 
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1 Cm 

(2021 - 594777 - 488377^) 

756 r 

-i-(131 - 90777 + 127377^)^ 7^2 
756 r 



1 



— (355+ 1906?7 - 33777^) 
252 



(24 - 72r;) 



42 c2 
rr f 1 



(26 - 20277 + Amf)v'^ 



1 Gm\ 
— (1085 - 405777 - 1463772) — j 



21 c2 
1 



(11 - 3377) 



(41 - 33777 + 733?7^)7;2 



126 



-(l- 577 + 5r,2)f2 
63 

1 Cm 

— (742 -33577 - 98577^) 

189 r 



(22) 



^ JI. STF 



ijk ' 



1 + ^((5 -1977)1-2 -(5 -1377) 



rr 

-(1-277) 



Gm\ 
r J 



-(I-277) 



(23) 



(1 - 377) 



110 C2 



(103 - 7357/ + 1395772)^2 



(100 - 6IO77 + IO5O772) 



55 c2 
78 r2 
55 c 



(1 - 577 + 577"^) 
2 (1 - 577 + 5772) 



/i— (l-2r,)STF,,fe,„{a;'^'='™} 



/i(l - 5?7 + 5772)STFijfc;„„ { 



ijklmn 



} 



(24) 

(25) 
(26) 
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The current-type moments read 
5„ 



1 



1 



28 c2 
-(54 + 60ry) — 



(13 - %%T])V^ 



28 c- 



:(5-10?7) 



J. 



ijkl 



ni / 



(27) 

(28) 
(29) 



The notation STFy... denotes the symmetric trace-free part of the tensor with indices ij... . The 
GW ampHtude, and from that computed, the far zone angular momentum and energy transport, 
has been completed to 3PN in [271 • From Thome's paper [T7], see his Equation (4.3), we 

also extract that the GW amplitude can equivalently be expressed in terms of tensor spherical 
harmonics, 



. TT 



1 



OO 

1=2 m=-l 



1^2 (I) 



r''{t-R/c}T: 



E2.lr 
3k 



l-l {I) 



S^"\t-R/c)T^, 



B2,lr: 



(30) 



The components l'™ and 5'™ are projected out of Equation P7)) with the help of the spherical 
components Y^™. The basis transformation and the explicit representation of the Y^™, taking the 
direction of the orbital angular momentum as the preferred one, h = /le^, reads 



S''"(i) 



167r 



(2/ + 1)!! 
-327r 



(Z + !)(/ + 2) 



1/2 



2(^-1)/ 

■(? + !)(/ + 2) 



1/2 



ii + i)i2i + iy.\ [ 2{i-i)i 

2/ A- 1 

(-1)"(2I - 1)!! 
(-1) 



ATi{l-m)\{l + m)\ 



3 



21 + 1(1- my. \ ' ^ (—1) (2/ 



2/c-l)!! 



Y. 



A, 



An {l + my.J ^ {I - m - 2ky.{2k)U 

(-1)™y1|"I* for (m<0). 



*2fc-f-Tn + l 



(31) 
(32) 



(33) 
(34) 



where the notation "[x]" denotes the integer number which is the largest smaller than or equal to 
X. From indices comprised by ( and ), the STF parts have to be taken. The number I tells us what 
type of moment we have: quadrupole for I ~ 2, octupole for I = 3, hexadecapole for Z = 4, . . . . 
From Equation (j34p we see that we only have to compute the moments for m = 0. . ./. Equal 
sign B above, taken from Appendix (A3) of |23) . may be more practical for programming than A. 



The relevant pure-spin tensor harmonics, T, 



E2,lm 
jk 



and T^^ 



are given in our Appendix [Cl The 
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reader can find the IPN outputs in [20], wliicli we also used for our prequel paper. In case of 2PN 
accurate conservative dynamics, the GW amplitude, consistently worked out to the same relative 
order in = e and having moved to our scaling, explicitly reads 



R 



nE2.2r, 



2 (2) 

m=-2 

" 3 (3) 

Eg3mrpB2,3m 

.m— — 3 

4 (4) 

EgimrpB2,Am 

.7Ti= — 4 

5 (5) 



B2,2ti 



E2Ar, 



B2,5m 



■ 2 (2 

.m=-2 

4 (4) 

m= — 4 

5 (5) 

m— — 5 

6 (6) 

rn— — 6 



3 (3) 
m— — 3 



5 m 



6m 



(35) 



To reach there, we first have to compute the l^^ time derivatives of the STF mass and current 
moments. We obtain them by means of the accelerations for a compact binary in harmonic 
coordinates to 2PN order, a = aN + aiPN + a2PN, (also taken from [22| . written in their units), 



aN 

aiPN 
a2PN 



Gm 

— 5-ni2 , 



r' 
Gm 



Gi 



-2(2 + 77)— + (l + 37/)i;2_Lr2 
r 2 



15 



(36) 
(37) 



^(12 + 2977)^^ + 77(3 - 477)t;* + ^77(1 - 377)7^' 

-77(3 - 477)7;2f2 - 1^(13 - 477)— - (2 + 2577 + 277^)— 7^2 
2 2 r r 



ni2 



77(15 + Ari)v'^ - (4 + 4I77 + 8r]^)—^ - 3t]{3 + 2T])r^ 



(38) 



The results are lengthy, easily reproducible and needed only temporarily; thus, they are not listed 
here. We employ the orbital parameterization, Equations ([T]) - (|16bp . to obtain the normal and 
the velocity vector ni2 and v in spherical coordinates, symbolically 



ni2 
dr 

" dt 



ni2 7', 

{cos((/)),sin((/)),0} , 



dr 

du 



du 



DM 



dM 



at 



(39) 
(40) 

(41) 



Equation (|¥T|l is to be computed with the help of the KE © . Again, it is not necessary to provide 
the velocity as functions of u and (j) because these terms are easy to be reproduced. With this 
input we can compute the spherical tensor components of Equation (|35p . The next section gives 
the results, using the 2PN accurate QKP in harmonic coordinates. 



4 The GW forms: tensor spherical components of K 

Defining 

A{u) := I — et cos u , 



TT 



(42) 
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we get 



(2) 

j20 ^ 



\E\e 



2{r,~26)Wl . 3(3r;-l) , 2(ry - 26) 



7A{u)^ UA{u) 



7A(it)2 
(^252 We, (5 -2,) 



^(,(11, + 16) - 5) (e? - 1)) + ^^^y^l^^^y (((1877 - 323,), 



3682)6^ 



-756 We, (2, -5) + 323if + 1651, + 1666 



21^(u 



(686 - 2,(31, -f 197))e? 



-h,(41, + 625) - 1610 



W? (,(1517, + 7549) - 8645) 



126A(w)4 



+ 



We^^ (4,(79, + 179) - 217) 
42A(u)5 



(43) 



(2) 
t21 



= 0, 



(44) 



(2) 
t22 



V 5 I ^(u) A{uy 



l^{15rj-82)Wl + ^,iet{6rj-79)sm{u)We, , 



3(1-3,) , 3 



S, - 210 



14(et - 1) 



+et(35(5 - 3,) + 3et(5et + 2)(3, - 1))) - j^et (t^^^ ((3 - 9v)et^ + , + 5) 
+2i ((9, - 3)6*2 ^ -^g^ _ 25) sin(M)) 

Wl (614,2 _ 37;^o, + 427) We, (,(4049, + 3901) - 10016) 



iet sin(M) 



63^(u)5 189A(m)4 

((,(377, - 575) + 2200)e2-, (3024 Wet +449, + 1933) + 7560 Wet - 1132) 

63We,A{u)^ 

^— -r f (,(11, + 16) - 5)e* + 2(,(67, + 398) + 59)6^ + 252 We, (5 - 2,) 
\A(u)'' V 

^ + iL(5-,(ll, + l6)) 



2lWlA{u) 



61,2 ^gg^ _ + —(,(11, + 16) - 5) + 



42 



A{u) 



1 



126We2 A(u) 



6(,(11, + 16) - 5)e^ + (,(2063, + 7867) + 5038)6? 



-3780 We, (2, - 5) - 2129,2 + 2621, - 11056j 

(4(169 - 64,), - 766)e? + 504 We, (2, - 5) + 373,2 _^ ^^^j^ _^ gyg 

21A(w)3 

We2 ((1247 - 12113,), + 28061) W^ (16,(143, - 335) + 199) 



378A{uy 



126A(ii) 



(45) 



Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



11 



(2) 
^20 



0, 



((675 - 3639?7)e? + 3695r;-73l) 4(7; - 1) 



28W2 



A{u) 



-A{u)^^^^Wlp2Ab - 209677) + We, (2054?/ - 2077)ef siii(u)} 



(46) 



(47) 



(2) 








0, 


(48) 


(3) 






j30 


= 0, 


(49) 



(3) 
/31 



-i<p 



-6- 



Qiet sin('u) 



1 



-\EV 



3771 I 5A(ii) A{uY 

((57/ - 7)e2 + 19?/ - 17) - i We, (57/ - 7)et sin(?i) 
10 We2A(u) 



(57?-79)e^- 457? + 119 -4 We, (7?- 14) - 2ii(47/ ~ 17)et sin(7/) 

20W2 A(u)2 + 3We,A(?7)3 

3 W^ (277 - 23) - 2i We, (77 + 10)et sin(it) 
5A(?/)4 



(50) 



(3) 
j32 



= 0, 



(51) 



(3) 
^33 



-207W||i?|^/^We,r/<5„,e-3'^ 



A(t 



8We, ( We, +^et sin(u)) 
5A(u)3 



6 _ 2iet sin(n) 
' 5 5 We, 



(79 - 577)e?+45r/-119 



20W2,A(i 



z We, (5?7 - 7)et sin(7x) - 3(5?/ - 7)e^ - 577/ + 51 
30W|,A(?/) 

- {45 We,A(?/)3}"^ |6 We, ((5r/ - 7)6? + 777/ - 205) 



+2 



iet sin(u) (3(577 - 7)6^ + 14977 - 208) | 



We, (We, (3477 -191) + 147(77 -14)etsin(u))' 




(52) 


15A(77)4 


= 16y^|i?p777We,77(377 -l)ei< 


3111(77) 

A(77)3 



-\E\e' 



'im{u) ((87/2 _ 4_ 4) g2 107/2 - 77/ + 2) 

3W2,(3?7- 1)A(77)3 



2 (27/2 _ 2427] + 79) sin(77) W^, (2?/^ - 24277 + 79) sin(77) 



(9 - 2777)^(77)^ 



(3 - 97/)A(7t)5 



(53) 
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(3) 
^31 

^32 





32 [2^ 

+ \E\e' 



et sin(w) 



sin(M) 



et ((8772 - IYtj + 4) + lOtj^ - Itj + 2) 
12W3 (3r/- l)A(w)3 



(2772 - 2427/ + 79) et We, (22772 + 338?7 - 109) et \ 
18Wet(37/- ^ (12 - 3677)A(w)5 J 

i (15(1 - 37?)2e2 + 5 (19 _ 3) ,^2 _ ^3 yi;2 + 7) _ yya ^ 25) 

10W2 (377- l)yl(ti)3 
i (197772 - 23977 + 55) i Wl (597/2 „ gg^ ^ 25) 



(547/ - 18)A(w)4 



(6 - 187/)A(ii)5 



(54) 



(55) 



(3) 
533 



0, 



(56) 



(4) 
j40 



8 /^|p|2 „ J 5 - 5e2 6 5 
21 V 5 Ml - 377) - - 



-|i?|62 



9 Wi (207/(77/ - 128) + 831) 3 W? (207/(77/ - 128) + 831) 



110(377- l)A{uf 



(6677 - 22)A(u)4 



>V2 (577(277 + 163) - 218) 57/(827/ - 1699) + 2711 
55(377 - 1)^(m)3 ^ 55(3?/- l).4(u)2 

6(577(877 -41) + 52) 6(577(877 - 41) + 52) 



55(37/- l)yl(M) 



55(37/- 1) 



(57) 



(4) 
t41 



= 0, 



(58) 



(4) 
7-42 



i^V2^|i.|2,7.77(377 -l)e-^|l-^ 
27 We.et sin(u) - 2e? + | 



>V2 + i>Vg,etsin(u) 
2A(w)3 



A(w)2 



+ \E\ 



{88(377 - l)^(u)5 } ^ We, { -22 Wt, (3O772 - 25O77 + 79) 
5i W2 (77/2 - 44777 + 151) et sm{u)^ 

W2 (9287/2 - 116207/ + 3743) + i We, (2147772 - 103317/ + 3175) et sin(i 

132(37/- l)A(7i)4 

-75 (53W2, -66) 772 



2907?2 - 23577 + 58 



+ (330W2,(377 - 1)A(7.)2)- 

+30 (17 Wl - 22) ?7 + 129 Wl - 330 
+12i We,et sin(u) ((4O772 - 20577 + 52) 



495(1 - 37/)2e^ + 8345 W2 7/^ 
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-14285 Wlrj + 3902 Wl - 4957f - 2310r; + 825^ 
+ {660Wli3r] - l)A{ufy^ (22(e2(l5 {29Wl " 20) 



-20 (75 Wl - 2) ?7 + 459 Wl + 20) 

+30(1 - 3r]fe'^ + (30 - 675 W'l) rf + 140 (13 Wj^ 539 - 50) 

+3i We.et sin(M) ((98577^ - 5715?? + 1781) - 3257^^ + 4835?^ - 1561)^ 

4077^ - 20577 + 52 4077^ - 205?; + 52 
55(377 - l)A{u) ^ 55 - 1657/ 



(59) 



(4) 

= 0, (60) 
3 V 7 ' ' ^ " ' \^ A{u) A{uY 

2 I 43 



2i Wet et sin(u) + 9e'^ - 9 4i We, et sm{u) - 8e4 + ^ 
2A{uf ^ A{uy^ 



e^\E\ 



sin(77) 



i We,^ (9237/2 - 1368377 + 4523) e* 



44A(u)5 

i We.et (48 (4O772 - 20577 + 52) + 5633577^ - 105 1 757/ + 29219) 

330A(w)4 

iet ((98577^ - 57157/ + 1781) - 3257/^ + 48357/ - 1561) 

4iet ((4O772 - 2057/ + 52) + 29077^ - 23577 + 58) \ 
55We,A{uy ) 
3Wl (7I6772 - 44567/ + 1461) 
UAjuf 

Wet^ (96 (4O772 - 205?7 + 52) + 14774077^ - 49688O77 + 149503) 

660A(w)4 

22 (et + 1) 



+ (330 We, A(m)3 {et - 1)) 'I -3e? (^9077^ 



/6072(et + l) \ 9900 (et + 1) 



^57/ ^ - 2363 We. j + - 3826 We, 



+1980 (377^ - 467/ + 15) y'^ Y~^4 + 3 We. (24307/^ - 18545?7 + 6074) 

+ We, (284 1 07/2 - 837957/ + 25262) et + We, (-284 1 07/^ + 8379577 - 25262) 



+ (330 We,A(7.)2 [et - 1))-^ I -6e?(5772 {^^^^^^ 
+77fl330We,-^^lMUl5%±l)_,9We, 



827 We, 



We, y We, 
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825 (et + 1) 



-416W, 



5e 



et ] -r IO77 



231 {et + 1) 

We, 



164 



V We 



(594772 + 277277 - 990) 



We, (-662877^ + 540777 - 1165) 



-2970(1 - 3?7)^ 



2 /l±£i, 



4- 5 We, (6O34772 - 81797/ 2155) 



-4O772 + 205?? - 52 1 



55yl(w) 



— (4077^ - 20577 + 52) 
55 



(61) 



(4) 
^40 



= 0, 



(62) 



(4) 
^41 



-y2^|S|5/2(l- 277)775,„e~*^ X 



W„^ 



2 We 



3A(u)4 5A{uy 



2iet [el — l) sin(u) 
3A(w)4 



(63) 



(4) 
gA2 



= 0, 



(64) 



(4) 
543 



= 40 



-|i?|^/2(l- 277)77<5„ 



W^ 



2 We 



■iA{uY 75A(w)2 



27et (gj — 1) siii(u) 
15A(w)4 



(65) 



(4) 
544 



= 0, 



(66) 



(5) 
j50 



= 0, 



(67) 



(5) 
j51 



2 /7^. 1^,5/2 ,^ '^l . 8^ We, 4»We, 

3 V n '^^'^ - ^) I(^ + TiR - 5W 



4iWe, , . , J 2et 
+ sin(u) 



3A(u) 



3A(w)3 7A(ti) 



(5) 
j52 



(69) 



(5) 
j53 



||i?|5/2,7(2,7-l)5„e-3'^x 
fet sin(u) - 8i We, 28i We, (e? - l) (16et sin(M) - 5i We, ) 



A{u) 



5A(m) 



3A(w)4 

2 (et (41 - 48e2) sin(u) - 6i (e* - 1) (8e| - 1 

9A(m)3 



(70) 



(5) 
^54 



= 0, 



(71) 



(5) 
j55 



27r 



= SW — |i;|5/2^(27?-l),5„e- 



■510 



165 
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16i ( We, + iet sin(u)) 3 (e^ - l) (16et sin(u) + 15i We,) 

i (-10 We, (24e2 - 17) + iet (l37 - 1446^) sin(M)) 
12v4(m)3 

7i We, -et sin(u) + 5i We, 1 



(72) 



(5) 
^50 



16 / 57r 

TV 231 '■^'^''^^'^^ " + ^ 
f sin('u) 2 Wet sin(u) 2 sin(u) 



A{u) 



5A{u) 



3A{u) 



(73) 



(5) 
^51 



0, 



(74) 



(5) 



f\l^\E\H5{rj-l)r^ + l)r^me-^^^x 



iWi 



iWl 



A{uf 3A(w)4 5A(u)3 



6«We^, , . / Wl 

+sin(u)etW, ' 



\A{uf bA{uY 2,A{uY 



(75) 



(5) 
553 



= 0, 



(76) 



(5) 
^54 



16 / In 

Y^Y^I^^I ??(5(??- l)r/+l)me 



—4:i4> 



35i Wt, 



6i {e\ - 1) hi (e? - l) 



1 A(uf 

+ sin(u) Cf Wet ~ 



^(u)3 
' 35 W^ 



3A(u)4 
1 



2A(uf A(uY 3A(u)4 



(77) 



(5) 
555 



0, 



(78) 



(6) 

j60 _ 



64 /IOtt 



33 V 273 

I- ' 



|£;|^?7(5(7/ - l)r/+ l)m x 



77 



IWt 



IWt 



A(u) 60A(m)2 10A(m)3 24A(m)4 8A(u)5 



(79) 



(6) 
^61 



0, 



(80) 



(6) 
j62 



'^'':|i?|V(5r]'-5r/+l)e-2.0 



33 V 13' 



1 |^-2e^ , 23 We^, 



W2 



107 Wf^ 

A{u) Aiuf- ~^ 24y4(u)5 ~ 10y4(u)3 ~ 72y4(u)4 

2 2 11 



+iet Wet sin(u) 



W2 



12A(m)5 A(m)2 3A(u)3 9A(ii)^ 



(6) 
j63 



(81) 



(82) 
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(6) 
t64 



(6) 



1 - 



^-^|i?|V(5r?'-57y + l) e 

1 - 8e2 (et - 1) {et + 1) (STBe? + 167) 



A{u) A{uY 72A{u] 



4 



2z>Ve,et (6et-5) (6et + 5) 35m^e/ 



(83) 



(6) 

/'^^ = 0, (84) 



27r 




3 13 yi;2 y^;2 (ll52e2 _ 551) 



8A(m)5 2A{uY 24A(w)^ 
We, (96e2-85) 2We, , 6 We, 



3A(u)4 ^(u)3 A(w)2 



(85) 



It will be necessary to decompose these tensor components in terms of irreducible expressions to 
get a time Fourier representation. Those will be terms which collect contributions having sin u on 
the one hand and those without sin it on the other, and they will be used when we write down 
the exponential of the orbital phase in such a way that we can use results and representations we 
already know from the literature or we have to evaluate them from scratch. This computation is 
subject of the next section. 



5 Relevant Kapteyn Series of irreducible components 

5.1 Series representation for the inverse KE, sinmw and cosmit 

We recall the computation of the IPN version of this consideration in paper 1. There we required 
only a Newtonian accurate expression of the sin- and cos-function of multiples of the eccentric 
anomaly. For further considerations, let us call u = g{M.) the solution to the 2PN KE. As 
we Taylor expand the argument of the Bessel integral - which will be done below - we are 
in the position to provide w as a series in Ai up to 2PN. Therefore we need the well-known 
representation of v and u in the KE ([2]), where in the 2PN term, we can insert their Newtonian 
accurate summation surrogates (see Equations (5) on p. 553 and (8) on p. 555 in [30] • In EI] there 
is a misprint in the definition of the Gn on page 33: the factor ^ should comprise the complete 
right hand side), 

oo 

W(e=e,) = M+^G,n{et) svamM , (86) 

m— 1 

G,„(e) = ^ -j J,n{me) + ^ a'* [J„i-s{me) - J,„+s(me)] |> , (87) 



=1 
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and a is extractable from 
2a 



1 + q;2 ■ 

The series expansions of the functions {v — u) and sini; at Newtonian accuracy read 



{v-u) = + ^G,(et) sin(i7W)j - + ^ -J„(net)sin(nA^) 

oo 

1 - ^2 J„{net) sin(n Al) . 



(89) 
(90) 
(91) 

We take above definitions and write in shorthand notation for further calculations, cf. Equation 



M = u — ct sinu + oij sin(jA^) . 



(92) 



Inserting this in the KE and solving for the Fourier-Bessel coefhcients, we calculate after Taylor 
expansion in e (see Appendix 



g{M)-M = ^Xsin(nA^), 

n=l 

2 r 

An / cos{nM)dg{M) 

mr Jo 



(93) 



2 

mr 

2 



cos n 



dg{M) 



9n{M)- et sm{gN{M)) + X! sm{Mm) 

m—1 J / 

cos {n [gNiM) - et sin(.gAr(A^))]) dg{M) 
-/ 8in{n [gN{M) - et sm{gNiM))]) e'^ amSm{Mm)dg{M) . (94) 

TT In — ' 

•^^> m=l 



Defining 
Q{j,n) := 
the result reads 
u = M + 2 



0, J<n 

1, j > n 



(95) 



E 



sin(jA^) Jj (jet) 



-2e^T,,^l- el ^ sin( ) 



Jra {etm) J'j+rn {et{j + m)) 



(etm) J'm_j (et(m - j)) 

.m=j+l 
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Jm (etm) J' -j-m (et(j - m)) 



.m— 1 



e(j, 1) + J', {ea) 



i-i 



.m— 1 



-Jrn {met) ( Gj-„j (et(j - to)) + 



2Jj-m (et(j - to)) 



("^e^) ( Gm-j (e*) 

^ J„ (TOCt) ( Gj+rn (et) - 



.m— 1 



J - TO 

^ 2J„i^j {et{m-j)) 
j - m 

j + m 



e^i) 



-G, (et) + 



2J, (jet) 



(96) 



The reader should keep in mind Appendix C and D of paper 1 where care is taken of products 
of infinite series of sin jA^ and cos jA^ with arbitrary coefRcients aj and This consideration 
is necessary to collect for terms with the same positive frequencies in the above expressions and 
many more. 

Now we proceed with the trigonometries of mw, m £ N. We know from the symmetry of sinmw 
and cosTOU, that only sinjAi or cosjAl can contribute. Thus, we decompose 



^cr^sinjAl , 



3 = i 



COS mu = 7™ cos jAi . 

3=0 

The coefficients ct™ and 7^' can be computed using 



7™ = — / sin mu sin jM dM . 



COS mu cos jM d/A 



(97) 
(98) 

(99) 
(100) 



Switching from the integration over dAi to du in the above equations using the 2PN accurate KE 
and Taylor expanding everything to e^, we can perform the integration. One technical - but easy 
to manage - issue is to re-convert the arguments of the integrals for a simple application of the 
Bessel integral formula, 



1 r 

Jy{x) = — / dMCOs(j/u — a;sinu) . 
7^ Jo 

Appendix iBl provides the calculation. The results read 
sin(TO "") = ^ 0-™ sin j M + e"^ "^af sin jM , 

n=l j=l 

-rn ._ ^ {J._,„(j et) + Jj+m{j Bt)} , 
J 



(101) 



(102) 
(103) 

Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



19 



4 

a'; 



-J 



^ oo 

J j-\-m — n (^f U + ")) 

(et(?^ - j)) + Jj+rn+n {et{j + "-))}, 



cos(m u) = cos j + e"' ^ 7^* cos jA^ . 



■= J yo-m{j et) - Jj+r-iij et)} X e(j, 0) + 5,nl Soj (-y ) , 

4 1 ^ . 

■= ^ Q!„| J-j+m-„ (et(-0' + n))) + Jj+rn-n {et{j - Tl)) 

n=l 

-J-]+m+n (etin - j)) - Jj+m+n {dtij + 'l))} X e(j,0) 

+<5ml <5oj 



(104) 
(105) 

(106) 



(107) 



This prescription is valid in both ADM and harmonic coordinates, where, of course, the appropriate 
values always have to be included. 

5.2 A(w)~" and j4(m)~" sinw as Fourier-Bessel series 
We keep expression (AlO) from [T], 



A(u)-" = 1 + 4"' + cos JU , 



smu 



Aiu) 



- ^ ys 



(") 



sm ju . 



(108) 
(109) 



In the above two Equations, bj"'' and S'j"'' are expansion coefficients for harmonics of u. From 
Equations (jl02p and JTUS]) we learn that the Fourier-Bessel representation of the above two series 
get 2PN corrections,!^ 



00 / 00 



A[u)- = 1 + + ^ ^™ + E E ~W 

m=l j"=l \fe=l / 



= Y^Af cosjTW, 



(110) 



00 / 00 



sin w j4(i 



= E E'^^rM--^-^ 

00 

= ^5f)sinjA^, 



fc=i 



r,(») 



(111) 

(112) 



"'^Please recognise the misprint in Equation (39) in paper 1: between feg"' and the round bracket, there should be a 
"+" instead of a 
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(") 



E 

k=l 



(113) 



5.3 Computing one more relevant double series 

We have seen in the IPN case that we have to deal with products of sine-sine and sine-cosine series. 
At 2PN accuracy, we additionally observe products of cosine-cosine series. The old double series 
formulas are valid irrespective of the PN order, but we supply the computation of cosine-cosine 
series products below. Suppose a term of the form 

/ oo \ / oo \ 



coskM E cosmTW 



\m—l 



oo oo 



k=l m=l 



Collecting for contributions with the same the frequency (for k — m ~ ±j we obtain m = k ^ j 
and for part 2, k + m = j we obtain m = j — k), the result reads 



oo oo 

- ^ ^ (cos [{k - 'm)M] + cos [{k 

k=l m=l 

= l(±A^k'B, 



i)M]) 



Kk=l 



cos jM X 



E A 

k=J+l 



(") 



k-j 



fe=i 



X:4"^s.-.e(j,i) 

fc=l 



CC,[n] 



COSJ.M. (115) 



j=0 



5.4 Decomposition of exp{— im0} 

What we have done at IPN accuracy has to be extended to 2PN, especially at the orbital dynamics. 
It is helpful to find a special decomposition of e"*™*^ in such a way that the mode decomposition 
of any IPN function of u (to be performed exactly) is not required at this point of calculation. We 
will combine the terms in such a way that we can use results known from the previous sections, 

{(p - (f)Q) = (1 + /C) V + (/40 sin2ii + g^^ sin3u) , (116) 











tan 








.2. 



V = v{u) = 2 arctan 



=: {M + wn) + viPN + e'' 'y2PN , 

-im((j> — <po) A g-'i m [(l+K;)ti + (e''/4^) sin 2-U+ (£^94^ ) sin 3u] 

A g-im u g-im(£^A;iPN+c*K:2PN) [■UN-l-e^i'lPN-(-«'''U2PN] g-i m c*" [/4,#. sin2ii+g40 sin 3-u] 
g-'imu g-imt)NC^('CiPN+'C2PNe^) g-i "i e"" [Kipn^ lPN-(-/4,/. sin2ii+g40 sin 3d] 
g-i m (e^K;iPN-|-e"*'C2PN) A-1 ^ 



(117) 



Part A {'in) 



Parts i'f^i) 
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[l - ime'*(/CiPNUiPN + /40sin(2uN) + 944, siii(3uN))] 



(118) 



Parte ("i) 

In step 3, vn decomposes into a purely secular part, namely M, and a purely oscillatory one with 
zero average over the interval Ai = [0, 27r], viz. wn [see Equation ([55|) ]. Therefore, in step 7, the 
middle term in edgy brackets ("Parte") can be written in terms of single and double summations 
of terms in the right hand side of Equation (|86p . almost without computational cost. PartA and 
Parte will contribute A{u)~" and also A{u)~" sinu terms that will multiply with the series of 
Parte. These contributions are lengthy and we again skip explicit expressions. In principle, other 
decompositions are valid as well, but we find it convenient to choose the above one because of its 
structural clearness. Parte is essential and we give it explicitly, 

1 — im/CiPNVNe^ — —mv^e^ (to/Cipn^Cn + 2z/C2Pn) 



Parte (to) 



= 1 — zTOe" K. 



IPN 



5Z '"^ipnGj (et) ■ 



fc=i L2 



to/CJpjj cos(fcA^) X 



^ G,(ei)G,_fc(et) 

j=k+i 



fc-i 



^G,(ei)Gfc-,(et) 



xe(fc,i) 



+4i/C2PN sin(fc7W)Gfe (e*) 



(119) 



J2J1 



For a clear understanding, we sometimes have added auxiliary indices to the brackets. This simply 
helps to see how deep the bracket in the current expression is. To IPN order we recognize what 
we computed in paper 1, Section IV. We will face products of Parte with powers of A[u) . They 
can be put into the form 

Parte (q) 



A{u) 



= ^ ^ige^ (/CiPN + IC2Pn€^) sm{jM) 

r j-i 



E 



-^(;^)e(j, 1)0,-™ (et) 



E 

m— j + l 



£ ~A^G,+m (et) 



24")g, (e,) + 



1 °° f 
--5]cos(j.M) ICjp^qh' 



^In Equation (C5) of paper 1, the very last term should get a factor ©(j, 1) to make it consistent with our notation 
and for convenience of the reader. 
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fc-l 



E 

k 



n=k+l 



fc-l 



E 

fc=i 

fc N 

-^[G„(et)Gfc+„(et)] 



n=k+l 



^e(j, i)E 



- E (^*) (^*)] + E (^*) (^*)] 



n— s+l 



^[G„(et)G„+, (et)] 



- [G,n (et) G,_™ (et)] + ^ [G„ (e^ G„_, (e*)] 

m— 1 n— j + 1 

j 



E^" 



.fc=i 



lElkiPN'^vf^^f- [G,(eOG,_,(eO] 
[Gfc (et) G,_fc (e,)] - £ [G^ (e*) G,+fe (e*)] j 



fc=i 



fe=i 



The part including sinw reads 
PartB(q) . ^ ■ / ■ ^ a\ 



■^0 



25, 



E 



1 / CO 

-/C We'e(j, E [Gfc(e*)Gfe_,„(ei) 



m — i m 

^ [-Gfe (et) G™_fc (et)] + [^^^ (e*) ^'^■+'" (^*) 



k=l 
oo 

E 

m— j-j-l 
m — 1 



1 / ^ 

-/CW^'^'^f^lJ Y [Gfc(e,)G,_,„(et)] 

\fc=m+l 



^ [-Gfc (et) G™_fe (et)] + ^ [g,. (et) Gfc+,„ (et) 



fc=i 



E 



2 „2^4<^(n) I 



(120) 



oo 

Ga; (e^) Gk^rn (et) 
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-Gfc {et)G,n-k (et) 



fe=i 



Gfc (ei) Gfc+m ( 



.fe=i 



OO -J 



(OO OO 



m— 1 



~2^^^^ (^IPN + /C2PN^ 



(121) 



This will be shortened by writing 
Parte (q) _ jsr^r a 



Parte (g) 



E(A)r''' cosjM + E(-^«)l"''' sin jAl 



E(5c)r'^' cos + E('5.)r sin^'-^ 

J=0 J=l 



(122) 



(123) 



with {Ac)^''''^\ (A)'"'''!, (5c)t"''?l and (^s)!"-"?! to be easily extracted from Equations ([MH) and 
and as well remembering that .4^"'' and ^j""* have 2PN terms. Equations (|120p and p2ip appear 
to be relatively complicated, but they are simply expanded applications of the product formulas 
for sin and cos series, Equations (CI) and (Dl) of paper 1, and our current Equation (|114p . As in 
paper 1, we can now decompose the wave form into the above irreducible components, from whose 
we have extracted the time Fourier series representation. To simplify matters, we introduce 



Parte (g) 
Parte (g) 



= : ^Slnq]{u) ■ 



(124) 
(125) 



5.5 Multipole moment decomposition: a brief posting of the results 
Using the results of the previous section, we decompose the multipole coefficients / and S as 



(a) 
jam 



Ct „--mHCM -miipo 
^1 -.am ^ c y\ 



Qtjra ri ^ — miKM „—mi<hQ 

J — L,S;bm e ex 



E'' ^Q;[fcm]-^[fcm](") + E'' ^^lkrn]J^S[k'm]iu)} , (126) 
. k k' J 

E''"'/3[fc„]-F[fe„,](u) + Y,''-^P[k^]Ts[k'^]{u) \ , (127) 



symbolically, where Ci-am and Cs-am are some pre-factors and a, a, (3 and /3 are coefficients to be 
determined, k and k' are some summation dummy indices with boundaries depending on a and m 
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and depending on the type (/ or S), and a — 2 and 5—1 are labels for a, a and for /3, /3, counting 
the order of c"^. The associated components of / and S, including PartA and Parte, £^re given in 
Appendix |D] Some "pre - Fourier" domain reads 



(a) 

= C/;ar„ e-'"^'^^ e-'"^'^" X <i > ' sin jVW /s" + > ' cos jTW /g"' *> , (128) 



(b) 



l^sinjAl /|f+^cosj.M /|H 

( (b) (b) ^ 

Cs;bm e-™'^^ e-™^« X 1^ sinjX 5^™ + ^ cos ^^^^j , (129) 



with 

}al (^-2^[,,„,(A)f"l+^"-^5[,,„,i('5.)f''"'V (130) 



^ k' J 



1% := I > '°~^«r.™,(A)r''"" + > ''^-^5rH„,i(5.)r I , (131) 

4t ■■= fE'^'AfcH(A)?'™l+E'"'%H('5.)f'"'y (132) 
\ fe' / 

■■= fE'"'%™i(-^4''"'+E'"'A^'™]('5c)?'"') , (133) 

\ k k' J 

for a e [2, 6] , 6 e [2, 5] , 



and for extracting the pure Fourier domain representation with delta distributions - (and not the 
one mixed in exponential and trigonometric representation as in Equations (jl28p and (|129p ) - we 
take the Fourier transformation of the sinjA^ and cosjA^ terms, 

r e-™'^^ sin jTW e*"*dt = 5 {jN + KmN - uj) - 5 [jH - KmN + uj) {IM) 

V27r j_oo V 2 V 2 

f e-™'^^ cos jTW e^'^Mt = ^5{jU + KmN - uj) + ^5{jU - KmN + uj) . {lib) 
This is valid only in the purely conservative orbital dynamics case. 



5.6 The effect of radiation reaction 

Up to now, we have used only the conservative parameterization of Section [2] hitherto. When we 
include the slow-in-time variation of the orbital elements due to radiation reaction, the Equations 
(|128p and p29p are no longer valid, because the radial frequency M is no longer a constant. In 
paper 1, we gave a detailled overview over what to do when eccentric GW signals with chirp have 
to be Fourier transformed using the algorithm of the steepest decent or equivalently called "the 
method of the stationary phase". The only difference to be recognized is that we may include 2PN 
accurate formulas for et and M and the QKP, but the analytic integration of those EOM will be 
skipped because of the extensive space filling of the results as we have faced already at IPN. A 
future publication may extend the calculation of Section VI of paper 1 with the help of Appell's 
integral formula (see e.g. [32]), the orbital averaged 2PN accurate EOM 

(136) 
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''^^') = e2PN(AA,eO, (137) 

orbit 



dt 



to be taken from the expressions (4.32) and (4.34) of reference and converted appropriately to 
Af and et, and the perturbation routines for inverting the formulas for the time of steepest decent 
we are familiar with. 

The interested reader may find the stationary phase approximation (SPA) results for circular 
orbits in various data analysis papers [331 1341 1351 136j . 



5.7 Considerations about limiting of the series 

Infinite series have to be restricted to finite ones for practical issues. In paper 1 we gave instructions 
how to limit the IPN series as we set up the following properties concerning the orders of involved 
terms: 

1. Bessel functions of order n: J„(nef) ~ ©(e"), 

2. The v(u) expansion coefficients: Gn{net) ^ 0(e"), 

3. Even-in-u expansion coefficients: A'j'^ - O(e^), 

4. Odd-in-u expansion coefficients: 5^"' ~ 0{e{~^), 

5. The sini; expansion: J^iet) = 5 (-'^ri-i(et) - J„+i(et)) = 0(e"~^), 

6. Double series expansions: 

These computations are still valid at 2PN and have to be applied to each series where ef"™ mdcx 
plays a role. The result of this is obvious but lengthy, so we skip the provision. The interested 
reader may take a look at |30| for more information on an estimate of the error when using 
finite sums. Having double, triple, . . . , maximally ri-tuples of summations, each evaluated up to 
some order 0{ef^) and, thus, containing DJl terms (plus or minus some finite number), we have a 
computational cost of ~ 9Jt" terms per each time step Aik and NSP x 2Jl" term computations in 
total, where NSP is the number of sampling points, k = (0, . . . , NSP — 1). In our case, a typical 
value could be 71 = 5. 



6 Conclusion 

In this article we provide 2PN accurate GW forms in terms of tensor spherical harmonics. The 
analytic Fourier-Bessel series of the inverted KE as well as the sines and cosines of the eccentric 
anomaly have been evaluated. These series may be interesting for perturbation theory of celestial 
mechanics. We could verify parts of the results of [5D] and reproduce results of the ideas of paper 1. 

One slight drawback should be mentioned. Without appropriate optimizing, the presented 
routine is far less than quickly computable, for example for a data rate of 4096 points per second 
and several minutes to be observed and a restriction to errors of, say, < 0{el). A first numerical 
insight, done in C for the IPN case, showed a CPU time of ~ one minute for the case of 128 data 
points at an error of 0(6^). It can, for example, give an impression of the orders of magnitude 
of how many harmonics may have to be included for a data analysis investigation. Its CPU time 
consumption should be improved by atomising the series computation to make it attractive for 
researchers in data analysis. 

A future investigation may include the 3PN GW amplitude ^\ and the 3PN QKP for point 
particles without spin, or even more the spin dependent multipole moments of |37l [551 [5^ HP] 
and |41| using the QKP for aligned spinning compact binaries and a well-suited and optimized 
numerical implementation. 



Copyright line will be provided by the publisher 



26 



M. Tessmer and G. Schafer: Full analytic frequency- domain GW forms 



Acknowledgements We thank Jan Sperrhake for useful discussions. MT likes to thank Johannes 
Hartung for an enlightening introduction on automatise and efficient programming of Mathematica, in 
which we coded most of our work. This work was funded by the Deutsche Forschungsgemeinschaft (DFG) 
through SFB/TR7 "Gravitationswellenastronomie" and the DLR through "LISA - Germany". 

A The solution to the 2PN accurate Kepler Equation: an alter- 
native way to derive it and a useful check 

The 2PN accurate KE, 

A4 = u — ct sinu + {j^v-uiv — u) + sinf) 

=: u — et sinu + J^4(m) , (138) 

can be inverted by defining 

g {M - e-^ Ti{u)) ^ u . (139) 

as the appropriate solution, however it will look like. We can Taylor expand it around the New- 
tonian solution, 

gN{M)-g'j^{M)e^:F^{u)) ^u. (140) 

The Newtonian solution g is known, 

oo 2 

.9n (n et) sinnA^ , (141) 

^ — ' 71 

)i=i 

and reads 

u{M) = gN{M) - ( 1 + X] '^Jn(net) cosnM] t^T^ {9n{M)) . (142) 



We know 



V — u = 

n=l 



E 



2 

Gn{net) Jn{net) 

n 



sin nM, (143) 



a/1 — e? sinw / ~ , , . , 

sinv = ^ = Jl-e^2y J'(net) sinnM. (144) 

1 - 64 cosu V * ^ ' ' 

n—l 

Inserting p43p and (|144p into (jl42p and applying the product rule for a sin and a cos series, we 
obtain the result of This is a nice calculation, so let us show it in detail. 
We abbreviate 

oo 

J^4{gN{M)) = ^ On sinnA^ , (145) 

n=l 

and read the coefficients a from Equations (jl43p and (jl44p . In fact, for the comparison, their 
form does not matter. Equation p42p together with the rule for products of a sin and a cos series 
reads 

E°° Jn{net) . 
sm nM 
n 
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-e^^ (a„ + 2P„^^'["l[j^,(fce,);afc]j sin nM , (146) 

{n—1 oo oo ^ 

e(n, 1) ^ a„_fc Jfe(fcet) + ^ a„+fe Jfe(fcet) - ^ afe_„ Jfe(fcet) I 
fc=l fc=l fc=n+l J 

n—1 oo 

= ^ afeJ„_fc((n - fc)ef) + ^ akJn-k{{n ~ k)et) 

k=l k=n+l 
oo 

fc=i 

just dropped out as for n = 1, for A: < vanish anyway. In summation with ak the last 
line gives ^^j^ ckfc ( J„_fc((n — A:)ei)— J„_)-fc((ri + A:)et)), remembering that for the k = n term. 
Jo(0)-l. 

Let us, in contrast, directly derive the expansion coefficients via integration and assume that 
g is the solution to the KE, u ~ g{.M)- Then, at mr [n £ Z), there are fixed points of the KE: 
mr = u ~ A4 and g{Ai) — M can be expressed in sin series, 



g{M)-M = ^^„sin(n7W). (148) 

The expansion coefficients, directly computed via integration read 

An = - I [q{M) - M\sin{nM)dM 
= — -I \g{M) - M]d{cosnM) 
= [boundary = 0] + — /" cos(n7W)[,9'(7W) - l]dA^ 



2 

mr 



cos{riA^} d{g{M)) 

— / cos |n — et sinit + e^J^4(it)) } dit 
'^7'' Jo 

2 2 ^ 

— / cos{ri (w — et sinw)} du / sin {n (u — et sinw)} \^ afe sin /cAl 

iTT Jo TT Jo 



mr 

2J„(net) 4I 



— / ak (cos [(k — n){u — et sin u)] — cos [(/c + n) (u — ej sin u)] ) du 
= ^^^^^-e'f]akiJk-ninet)-Jk+n{net)), (149) 



fe=l 



and we see that the 2PN coefficient shows agreement in both calculations. In step 5 we used 
Equation p38p , in step 6 we Taylor expanded the argument of cos around the Newtonian A4 and 
in step 7 we used that only Newtonian A4 is required in the sum. The trigonometries of u are 
dealt with equivalently. 
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B Fourier representation: sinmu and cos mu at 2PN 

As an exemplary calculation, we determine the expansion coefficients of cos mw, 

oo 

COS mu = ^ 7™ cos jM . (150) 

j=o 

Using integration by parts, the computation turns out to be 
2 



7™ = — / cos mu cos jA4 dA4 

= [cos mu sin jA^] J — — j (dcosmw) sinjTM dTM 

2m r N • ■ , ^ du 

= ... n / (sin mu) siniM —t-iOlM 

3^ Jo ' ' AM 

= — / (cos[tou + j7M2PN] ~ cos[mu — jA^9PN]) du 

= — / (cos[tou + jA^n] — cos[mu — jA^n]) du 

m °° 
/ (sin[mM + jA^n] — sin[mit — jA^n]) 7 sin /A^n 

^ Jo fr; ' 



/=i 



— {Jj+,nUet) - Jj-m{iet)) 



m °° 

— / a„ cos(jA^ — mu — Mn) — cos[jM. — mu + Mn) 

^ n=l ■'^ 



— cos(jAl + mu — Airi) + cos(jAl + mu + Ain) 
y {Jj+m{jet) - Jj-m{jet)) 

r2TT 



du 



{ 

n=l 



cos [(j — ri) (u — et sin(M)) — mu] 

— cos [{j + n) (u — et sin(u)) — mu] 

— cos [{j — n) (u — et sin(u)) + mu] 

+ cos [{j + n) (u — et sin(u)) + mu] |du . (151) 
The task is now to bring these integrals to the form 



1 r 

- / cos [x; (u - y; sin u)] = J^. (x^ y^) , 

7^ Jo 



(152) 



with some prefactor x and "eccentricity" y to be determined for each special case. In order of 
appearance in the last equation above, these eccentricities read 

yi =7^. (153) 

y2 =i^, (154) 

ys =fqM, (155) 
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(156) 



^ etU+n) 
^ j-\-m-\-n ' 

X,; is simply the denominator. Herewith, the rest is easy to calculate and so are the coefficients of 
sin mu. 

C The relevant tensor spherical harmonics 

We take the definitions from [T^, Equations (Al) - (A5) therein, 



LM 



LM 



A 



LM 



= Alm 
= 2Cl 



B 



LM 

Cl 

iE2 
LM 

iB2 
LM 



-2cJ|-cot 



- i Blm 

~ I Alm [0$ + 4,9^ , 



LM[ 



[2L(L + l)(L + 2)(L-l)]-i/' , 



(m) 



LM ' 



LM ' 



(157) 
(158) 

(159) 

(160) 

(161) 
(162) 
(163) 



with 6 and (jj being basis unit vectors in 9 and (j) direction {9 is the angle between the orbital 
angular momentum and the line-of-sight vector N, and (j) measures the angle from the x axis to 
N projected onto the (x, y) plane, see Figure 3 of [IB])- The relevant multipoles read 



E2 



E2 



^21 



rpE2 

^20 



rpB2 

^22 



rpB2 

^21 



rpLi2 ^ / 

on — ^ A / 



1 rs" 



,2ict> 



27r 

1 /15 . 2, 



-(cos(26i) + 3) + 2icos(6i) 



>s{9) sm{9) (§9 ~ + i sin(6l) (j)^ + (i>9 



-J-sm\9) 



^20 



T, 



E2 
33 

rpE2 

^32 

rpE2 

^31 

E2 
30 



-1 

Te V 27r 

1 /T 

4 V 27r 
1 /l5 

8 V TT 

1 [21 
"32V V 

1 rr 
32V 2^ 

1 /35 



20 



2i(p 



iff) 



32 
33 



32 V TT 
8 V TT V 

1 [21 



3id 



32 V TT 



(j)4> + 4>9^ (cos(26i) + 3) - Ai(j)9 - 00^ cos(6') 
(J)9 - +i(j)(i> + $9^ cos{9) sm{9) , 
sm'^{9), 

Ai(9<j) + ^9j cos{9) + [99 - <^(^) (cos(26l) + 3) sin(6l) , 
8i {94) + 4)9^ 008(26*) + (^^^ - (5 cos(6i) + 3 cos(36i)) 
Ai [94 + cos(e') + {99 - 4>4j (3 cos(26') + 1) sin(6') , 
-0^) cos(6') sin2(6l), 
(94> + 4>9j (cos(26l) + 3) - 4i - 008(6*) sin(6l) , 



(164) 
(165) 
(166) 
(167) 
(168) 
(169) 
(170) 
(171) 
(172) 
(173) 



(174) 
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rpB2 _ ^ I ■ 



-1 / 7 
32 V 2^ 
1 /35 



32 



,2i0 



.i<t> 



(5 cos(6i) + 3 cos(36i)) - 8i 6*6' - cos(26') 



B2 
30 



(3 008(26*) + 1) -Ai 
cos{e) sm^{e) , 



003(6*) sin(0) , 



(175) 
(176) 
(177) 



E2 
44 



7 



_E2 
43 



£2 
42 



32 V 27r 

32 V TT 



64V27r 



_E2 



E2 



T, 



rpB2 _ " I ■ 

J- A A — TTTT \ / 



B2 
43 



32 

32 V TT 

3 /JZ 
32 V 2^ 
3 [7 



+ cos{0) + [ee - (cos(26i) + 3) 8114^(61) , 

A[ee - (f>^^ cos^{9) + i (§(}> + (1)6^3 cos{2e) + 1) sin(6') , 
2i(e4> + (t>ej (003(6*) + 7008(36*)) 
(4co8(26i) + 7co8(46') + 5)j , 

(^0 + 0^)(7co8(26i) + 5)+ (^^^-0(/))(5co8(6') + 7co8(36*)) 8in(6l) , 
(7 cos(26i) + 5)81112(6*), 



32 V tt"" 



3i(f> 



B2 



T, 



rpB2 



(eij) + 4>0^ (008(26*) + 3) - Ai(j)e - cos(6*) 8111^(6*) , 

4(^9$ + i)ej cos^{e) -i{ee- ipif^ (3co8(26i) + 1) 8iii(6') , 

{e^j) + (4 cos(26*) + 7 cos(46i) + 5) 



64V27r 

-2i [ee - (t)<j) ] (008(6*) + 7 co8(36*)) 

-3 



T, 



B2 
40 



32V^ 
-3 fb_ 
32 V tI^ 



6*0 + (f)ej (5 co8(6i) + 7 co8(36i)) - i 
+ (^0)(7co8(26i) + 5)8in2(6i), 



(7co8(26*) + 5) 8in(6*) 



(178) 
(179) 

(180) 
(181) 

(182) 

(183) 

(184) 

(185) 
(186) 

(187) 



rpE2 

^55 



rpE2 

^54 



1 /165 , 



64 V TT 

64 V 27r ' 



E2 
53 



1 /33 



rpE2 

^52 



rpE2 

^51 



■5*-^ \ Ai{eij) + ij)e^ cos{e) 
i{e4) + 4>e^ (2cos(26*) + 1) 

(19 co8(6l) + 5 cos(36l))j sin2(6') , 

„ /''^[4if^0 + 00"V7co8(6*) + 9co8(36*)) 
256 V TT L V / 

+ [ee - (28 co8(26*) + 15 co8(46*) + 21) 8in(6*) , 
128 V 27r 

+ {ee - (14 co8(6*) + 3(cos(36i) + 5 co8(56*))) 

1 Pn_ 

128V 2^' 



(co8(26i) + 3) 8in3(6l), 



8i 6*0 + 06* (co8(26i) + 3 cos(46*)) 



(188) 



(189) 



(190) 



(191) 



.i<t> 



Ai [e$ + (5 co8(6*) + 3 cos(36*)) 
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E2 
50 



B2 
55 



B2 

54 



rpB2 

-^53 



32 
52 



-^51 



52 
50 



1 /1155 



(12 008(26*) + 15 cos(46') + 5) sin(6') , 
m - ^<f>) [5 cos(6') + 3 cos(36')] sm^{6) , 



[e^ + (cos(26i) + 2,)-Ai{ee- cos(6') siii^ (6*) , 
{eij) + 0^) (19 cos(6') + 5 cos(36i)) 



64 V TT 

^ /165 

64 V TT 

1 /33" 
^eiV 2^' 

-8i(^^- 0(^)(2cos(26i) + 1) sin^(6l), 

^ y^e^'"^ [(^^-^ + 0^) (28 cos(20) + 15 cos(40) + 21) 
-4i [ee - (^0) (7 cos(6i) + 9 cos(36i)) sin(6i) , 
1 [TT 



,2i4> 



128 V 2tt 
~8i(j)9 ~ (j)(t>] {cos{2e) + 3 cosine)) 

1 ffj 
"l28 V 2^" 



+ (jjO (14 cos(6i) + 3(cos(36i) + 5 cos(56i))) 



-">'^{eif) + 4>9^ (12 008(26*) + 15 cos(46') + 5) 
Ai (90 - (t>4>j (5 cos(6i) + 3 cos(36i)) sin(6') , 

— \ id<l> + (l)d] [5cos(6i) +3cos(36')]sin2(6'), 

64 V TT V / 



(192) 
(193) 
(194) 

(195) 

(196) 

(197) 

(198) 
(199) 



E2 
66 



E2 
65 



rpE2 

-'64 



E2 
63 



rpE2 

^62 



E2 



E2 
60 




256 V 27r 



4:i{9(j) + (j}0] cos{9) + [99- (j)(l)]{cos{29) + 3) sm^{9) 



- - (13 cos(6l) + 3 cos(36l)) sin^(6») , 
1 /l95 



1024 V TT 



8i 



(§4) + 4)9^ (13 cos(6') + 11 008(36*)) 



(92cos(26i) + 33cos(46i) + 67) sm'^{9) 



'13 



,3i4> 



2i 



+ (j)9 (52 cos(26l) + 55 008(46*) + 21) 



2id 



4i ( ^0 + 061 ) (10 cos(6') + 81 cos(36i) + 165 cos(56i)) 



1024 V 27r 

+ {99 - (126cos(6i) + 75 008(36*) + 55oos(56i)) sin(6i) , 
1 

4096 V 2^ 

i - 00 ) (289 oos(26i) + 30(oos(46i) + 7) + 495 oos(66i)) 
1 /65 

1024 V V 

+ (^^- 00)(7Ooos(6l) +87oos(36l) + 99oos(56l)) sin(6i) 



J4> 



2i 



(§4> + ij}9j (60 008(26*) + 33 008(46*) + 35) 



1365 



1 



512 V 27r 



[60 oos(26l) + 33 008(46*) + 35] sin^ (6*) . 



(200) 

(201) 

(202) 

(203) 

(204) 

(205) 
(206) 
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D The multipole components in terms of Fourier coefficients 



(2) 
j22 



>^[02] 



*[12] 



a[22] 



"[32] 



a [42] 



"[52] 



"[62] 



"[72] 



°"[12] 



X! °"[fc2]-7^[fe2] {u) + X! ""[fe2]-^S[fc2] (u) > , 



(207) 



. fc=0 



fc=i 



[2 - e?] + £2 -\E\ [(977 - 3)e? + Mrj ~ 442] + 



42 We 



(77(1177 + 16) -5)e^ 



+3(77(1977; - 1216) + 3717)e? + 32(40 - 63 WeJv + 5040 We^ - 6027;^ + 2798 



Ft - 2] +^'^^\E\ ((3 - 977)e? - 9477 + 778) + 



(5- 7^(117; +16))e^ 



+277(1512 We, + 3OI77 - 64) - 3e? (77(168 We, + 1977/ - 3184) - 420 We, + 10029) 



-7560 We, +2530 



(208) 



-2( We'j + ^\E\ (-977 (e? + 17) - 151e? + 163) + 



126 We^ 



€((11340 We 



-45322 - 77(4536 We, + 11237; + 15623)) + (3(951 - 119?7)?7 + 4188)e^ 



+27/(4536 We, + 74O77 - 2435) - 22680 We, - 57574 
1 



(209) 



:= 2[e?-l] +.^|i?| 



21 



Wl (11(377 + 13)e2 + 129?7 - 1163) 



\E\' 



189 



77(31752 We 



-343377 + 13700) - 79380 We, - 60259 + e2(?7(-4536 We, - 348177 + 10081) 
-11340 We, - 93512) + 3((830 - 4777)77 + 2117)e^ 



'\E\ 



21' 



WeVl577 + 121) 



378 



9(77(1177 + 3339) - 9453)e 



-477(6804 We, + 12117/ + 7423) - 68040 We, + 77866 



(210) 



(211) 



e'\E\ 



378 



3(47/(3497/ + 917) - 1117)e 



-1677(3087/ + 1891) + 95270 



Wl (7/(68537/ + 20333) - 12793) , 
-e^^l^l" We,'(477(5377 + 103) - 131) , 



(212) 
(213) 



i2We,et + i 



\E\e^et ((109 - 1977)e? + 477/ - 22l) 
TWi 



+ie" 



\E\' 



21 W|, 



-252 We^^(27/ - 5)et (e? - 2) + Wf^e* ((-4467/2 + 21447/- 4876) 
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°a[22] 



°a[32] 



°a[52] 



°<5[62] 



°Q;[72] 



+ {61{r] - 4)?7 + 1445)e^ + 77(301r; - 640) - 1399 



(214) 



l2\E\ef 



2 r6et ((54r; - 151)e^ + (600 - 14877)6^ - 4(8r; + 37)) 



7Wl 



~12(5-2?7)et 



-i2 Wlet - le^ —\E\ W^.e* ((15r; + 149)e? + 153?7 - 667) 



63 W,, 



3024 We, (277 - 5) + 7407;^ - 707?/ - 20795 



-64377^ + 165177 - 14042)e? + 5(77(317/ - 112) + 443)e: 



(215) 

-756We,(2r/-5) 

(216) 



IK 

189 



We.et 



6804 We, (2?/ - 5) + 12137/2 + 158217/ 



-2126 + ((8627 - 701?/)?/ - 37414)e^ 



^-\E\Wle^3rJ-l)et + ^e' 



+137/(317/ - 25) + 12269 



\El 
63 



We.et {el - 1) 3(477(1477 + 67) - 77)e; 



^\El 
189 



Wl^ {ri{3673r] + 14153) - 10828)et , 



ie^ ^\Ef Wl (477(5377 + 103) - 131)et 



(217) 

(218) 

(219) 
(220) 



(2) 



= 0. 



(221) 



(2) 
t20 



^[00] 



*[10] 



"[20] 



"[30] 



a [40] 



*[50] 



' fc=0 

1 - ^\E\e\3v + ^\E\' {-Ihf 16?/ + 5) , 



-l + -\E\e'(3v-l) + e''\E\' 



30 - 1277 1 , , 

+ — (77 II77 + 16 - 5) 

We, 42^'^ ' ' ' 



(222) 

(223) 
(224) 



-|£;|£2(77 - 26)+ , 
7I I V/ ) 126 W? 



-165177 + 1666 



((1877 - 323?7)?7 - 3682)e? + 756 We, (27/ - 5) + 32377^ 



-\E\e\2Q~ri)Wi 
-1610] , 

1 '-|2i.,2 



2 ' |£;p [(686 -277(3177 +197))e? + 7/(417/ + 625) 



-e- We, (7/(15 1 77/ + 75 49) - 86 45) 

We', (4^7(7977 + 179) - 217) , 



(225) 

(226) 
(227) 

(228) 



(2) 
^22 



= 0. 



(229) 



(2) 
^21 



'^\/T^^e^r]\E\^^'^ {mi - 7712) 
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{E ^Pm^mi^) + E ^Pm^smiu)] , (230) 

fc=2 fc=3 

|i^|e^((899 - 369577)e^ + 375l77- 955) 

A21] ^^y^l 1, (231) 

V[3i] := 

+ |£:|e2 ((675 - 363977)e? + 3695?? - 1067) + ^ (2077 - 2054?7)^ , (232) 

V[4i] ^|i;;|W,V'(103407y- 10497), (233) 

V[5i] := ^|i?|WeV' (2133 - 206877), (234) 

,~ i|£:|e2et((3639r;- 675)e2- 3751?7 + 955) , , 

i/3[3i] := -zWe,e, + ^ SsW^T ' ^^^^^ 

'/9[4i] := -i^|S|We,e'(2068r/- 2049)et, (236) 

'Asi] := ^i|i?|W,V'(206877- 2133)et, (237) 

(2) 

S-^^ = 0, (238) 



(3) 



/27r 



E ^"[fc3]-i^[fc3] (") + E ^"[fe3]-^S[fc3] (?i) \ , (239) 



fc=0 fc=l 

a[03] := -2et +8 — ^ , (240) 



1 _ ,2,o li^k" ((5r? - 7)e^ + (1 91r? - 8 29)e,^ - 26877 + 1124) 

6W2 

.2 ((7 - 57?)e4 + (1477 - 191r;)e^ + 4(677? - 497)) 



a 



[13] 



\E\e'^ ((237/ - 9Gl)e4 + (6305 - 220377)6^ + 254077 - 5764) 

36 Wl 

^a[33] := -e?(e2 + 3)+4 

+ ^\E\e^ (9(237; - 13)ef + (5353 - 39577)6^ + 5487? - 14668) , (243) 

ia[43] := 20W,4^ ^^l^lW^e' ((211^ +1162)e? + 3(3837;- 745)) , (244) 

ia[53] := -12 Wf^-i|i?|W4e' (9(297? + 38)e? + 3317;- 4973) , (245) 

'"[63] = -^|i?|W,V(14077 + 467), (246) 
90 

'«[73] = y|i?|W« 6^(477 + 1), (247) 
2iet(3e?-4) ijE^le^et (3(197? - 89)e? - 2687; + 1124) 

"['^1 ^= — w:^ ' ^'^^^ 



1 



1 ;r 



a[23] 



36i(77 - 4)et (3ef - 76? + 4) 36iet ((37? - 13)6^ - 4(77 - 5)) 



(249) 
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"^a[33] 



"^"[43] 



^a[53] 



"^"[63] 



^"[73] 



^ iet (29ef - 65e^ + 36) 



ef 



30i{r] - A)et {3ef - 7e^ + 4) 



iet (9e2 ((357? + 23)6^ - 2477? + 653) + 2 1407? - 6836) 



36 Wp 



+167 Wlet ~ 87 We.et {2et - + 3) + \E\^ [-7 (5^ _ 7)g^ 
+27 Wf^et ((577- 7)e2 + 1397; - 497) + h.W,^et{2{h^ - l)e\ 

+37(117/ - 40)e? - 5057/ + 1343)1 , 
12iWf^et 

+ \E\(^ (a&i ysHX iv " 4)et + ^7 Wlet ((1917/ - 472)6? + 737? - 479)^ 



16, 



i\E\Wle^i20r]+113)et. 
20 

— 7|i?|WX £2(477 + 1)6*, 



(250) 



(251) 

(252) 

(253) 
(254) 



(3) 
j32 



(255) 



(3) 
^31 



6 —6 "''"S\I'^\E\^''^ {mi-m2)vyVe,y- 



"[01] 



a[ii] 



^21] 



"[31] 



"[41] 



"[51] 



^"[11] 



^"[21] 



^"[31] 



X] ^"[fcl]-^[fel] (it) + X! ^"[fcl]-^S[fcl] (w) f" , 
1 

12' 



^-^AE\e^-'-^^ 

10 I I I yy2 



5?7-7 



j|£;|e2 ((5?/- 7)6? + 197/- 161) 



12 W2 

5i 5i|i;|e2((57/- 79)e2 + 197/ + 55) 
6" 



72 W2 



1 

24 



i|£:|e2 ((137/ + l)e? - 53?/ + 279) , 



6 

^ 7|£:|W2 6^(407/ - 317), 



18 



-7|i?|W,V (8^-49), 

et , |S| We, 6^(89 - 1977)6t 



We, 12 



6|g|>Ve,e'et 

hWe.et \E\e^et (3(79 - 57?)6? + 957/ - 157) 
6 ^ 72 We, 



(256) 

(257) 

(258) 

(259) 

(260) 
(261) 
(262) 
(263) 

(264) 

(265) 



Copyright line will be provided by the publisher 



36 



M. Tessmer and G. Schafer: Full analytic frequency- domain GW forms 



'5[4i] i|£;|We,e'(85-8r;)et, (266) 

i5[5i] i|£;|e2(8r?-49)etW3 , (267) 



P[32] 



8 /27r 



(3) 

3 V 7 

r 7 7 



E 'l3mJ'[k2] {u) + J2 'Plk2]J's[k2] (u) \ , (268) 



.fe=2 fe=3 



2i\E\e^ {{7]{lQi] - 61) + 20)e2 + (163 - 4677)77 - 50) 



%2, :^ — — 7^^(3:7-1) (2^^) 



2 , 2i\E\e^ 



2 



3(577(1677 - 25) + 28)e^ 



+2(5r/(7377 - 13) + 4)e^ - 577(19477 + 223) + 448 , (270) 

- 658)e2 + 
45(377-1) 



, 27|i;|e^ (3(577(7177 + 370) - 658)6^ + 577(66277 - 2477) + 3809) 
/3[42] := 107W,^ + ../o^„n ' (271) 



2i|£;| ^2^2 (-3 (-505^2 ^ 9^077 - 344) e2 - 577(2307; + 4147) + 6847) 



(272) 



2»|i;|>V4 6^(477(24577 + 604) - 845) 
A62] := , (273) 

27|i;|>V6 6^(2077(777 + 8) - 59) 
A72] := W—:; , (274) 



We, 

l^le^et ((77(877 -17) + 4)6^ + 3(77(2277 -95) + 30)6^-92772 + 32677 -100) ^ ^ 

, (275) 



(276) 



(277) 



3(377 - 1) W3 

'A42] := 4We,et 

2|i;|e2et ((10(38 - 577)7/ - 157)6^ + I66O772 - 1330?/ + 302) 
^ 45(377- 1) We, ' 

2;3[52] := -6etW3 

\E\e^{-et We,) (3(17077(77 + 2) - 131)6? + 11077(1077 - 67) + 2338) 
^ 15(37/- 1) ' 

2~ 8|i;|e2(477(3577 + 121)- 167)6t>V,^, 

/3[62] := , (278) 

2 2j 2|i;|62(2077(777 + 8)- 59)6tW5 

^= 977-3 ' ^ 

(3) r-i- 5 

= -16^ — |i?|'6,yr^/i(l- 377)^%o]-F5[fcol(u), (280) 

2« , \E\e^ (((17 - 877)77 - 4)6| + (7 - 1077)77 - 2) 
^[^"1 ■= ^ 3 W2 (377 - 1) ' 
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2|g|e^(2(?7- 121)?7 + 79) 

9(37?-!) ' 
\E\Wle\2{r^- 121)7^ + 79) 
3-97; 



(282) 
(283) 



(4) 
r44 



>^[04] 



[14] 



*[24] 



^[34] 



*[54] 



^[64] 



*[74] 



*[84] 



*[94] 



9 



9 



X] ^"[M]-7^[fc4](lt) + X] ^<5[fe4]-^5'[M('") f ' 
. fe=0 k=l J 



55(377-1) 



- [(5(41 - 87,)7, - 52)e 



-8(57/(12477 - 531) + 828)6^2 + 8(577(25677 - 1103) + 1708) 



(577(877 - 41) + 52)e: 



' ' 55(3?7-l) 
+8(577(12477 - 927) + 1488)e? - 8(577(25677 - 1895) + 3028) 

1 . . o ... \Ek^ 

6 



-29ef + 112e? 



-5(77(31477- 3185) + 1121)e^ 



' 330(377-1) 

20(?7(382977 - 12717) -f 3903)e? - 8(577(292177 - 8542) 4- 12788) 



- (3e^ - 19e^ + 24e? 



330(377- 1) 



(76IO772 - 7322577 -t- 22598) 



-3(1577(5877 - 69) + 218)e^ + 40(7?(411?7 + 10226) - 3407)e? 
-8(577(536774- 11245)- 18912) 



6) 



\E\Wle^ 



(2077(21477 - 2359) + 15563)et 



660(377- 1) 

20(77(970577 - 64529) + 20707)e2 4- 8(577(129537/ - 54357) + 84244) 



(607?(82477 4- 67) - 7371)e^ 



4-220(77(10777 - 363) 4- 145)e2 - 8(57/(756177 4- 45019) - 77612) 



140 yv° 



\E\Wle' 



3 9977 - 33 

-4(77(994077 - 26869) 8016) 



(377(451577 4- 11017) - 12361)e? 



lElW^e^ 



3(1577(21377 95) - 703)e 



-4(77(5357/4-34698)- 11300) 



6|£;| yVi°e2(77(122577 + 2427) - 844) 
3377-11 
70|£:| Wi2e2(,y(3577 -F9) -4) 
3377-11 ' 



(284) 



(285) 



(286) 



(287) 



(288) 

(289) 

(290) 

(291) 

(292) 
(293) 
(294) 
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^a[i4] 



^Q;[24] 



^'5[34] 



4^et {et - iej + 2) Ai\E\e'et 



^(5 [44] 



^"[54] 



^"[64] 



^5(74] 



^Q;[84] 



^5(94] 

(4r 

^42 



55(3?7- 1) We 



(5?7(5877 - 245) + 388)£ 



-(5(2165 - 50477)77 - 3364)e? + 10r;(256r7 - 1103) -f 3416 



24i|£;| We.e^et (e^ - Be^ g 

iet (9e6 - 55e^ OOe^ ^ 44) 
3Wl 



(295) 
(296) 



i\E\e^et 



330(377 - 1) We 



15(?7(133?7 - 815) 273)e? 



-55(77(100977 - 4619) 4- 1457)e^ 4- 4(577(828677 - 37381) + 58391)e 
-8(1157^(12777 - 578) -t- 20708)1 , 
i\E\We,e^et 



-7etWf^(5e?-12) + 



20 (42667^2 - 345277 + 829) e 



3 - • - / 330(377-1) 

-(1655(77 - 1)7; - 417)e^ - 8(57^(34577; - 2049) + 1796) 



-276* Wf^ (e? 



i\E\e'etWl 
220(377-1) 



(577(581?7 - 1637) 2473)e^ 



-220(277(7377 - 520) + 333)e2 -|- 88 (91577^ - 761577 + 2444) 

'ii\E\e^et (11(677(1577 + 37) - 89)e? + 4(5?7(34077 - 131) + 76)) 
^ 9977 - 33 ' 

-20zet 

i\E\e^et Wj^ (3(107^(19777 + 111) - 563)6^ + 4(7^(31407; - 24249) + 7712)) 

9977 - 33 

4i|£;|etWf^ £2(^(122577 -^3483) - 1196) 
3377 - 11 ' 
70i\E\et WeVe2(77(3577 -f 9) - 4) 
337^-11 ' 

\ ^"[fc2]-^[fc2] (w) + ^"[fc2]-^S[A;2] (u) > , 



. fc=0 



fc=l 



"[02] 



*[12] 



a [22] 



"[32] 



\E\e^ ((5(41 - 8?7)?7 - 52)e? - 2(577(12477 - 531) -f- 828)) 
55(377-1) ' 
2 lEle"^ ((577(87; - 41) + 52)e2 + 1077(12477 - 927) + 2976) 
^* ^ 55(377 - 1) ' 

lle^ \E\e^ ((632 - 5?7(17?7 + 478))e^ + 577(204777 - 3411) + 4797) 



6 



-3e^-t-lle?- 



\E\e^ 



330(3?7- 1) 
165(1 - 477)7/6^ 



g V . — t . 330(377- 1) 

-^(577(11867; - 9227) + 14738)e2 + 10(7369 - 52777)77 - 23978 , 



(297) 
(298) 
(299) 

(300) 

(301) 
(302) 
(303) 

(304) 

(305) 
(306) 
(307) 

(308) 
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a [42] .— 



2 ;r 



3 

\E\ W^^e^ ((875r/2 - 6815?/ + 2738) - 5?y(5715ry + 3613) + 8167) 



"[52] 2 KVg 



660(377- 1) 
33(r;(35r;- 191) + 60)e? 



55(377 - 1) 
6677 - 22" 



(309) 



^* 132(377 - 1) 

+(28933 - 622977)77 - 8915] , (310) 

|g|>V6£^(r;(437577- 29191) + 8993) 

"[^^1 ■= 132(377 - 1) ' 

|i;|>V« 6^(577(12577 - 653) + 983) 

"[-I ^= 44(3^FT^ ' 

^ 2i\E\e\t {{5rj{58rj - 245) + 388)e^ + 5(531 - 124r;)r; - 828) 
a[i2] := -2.W,e, + 55(377 - 1) We. ' ^^'^^ 

m\E\e\t (e? - 2) 

"[22] - yy^^ > (314) 

«[321 := + + 660(3,- 1) We. [15^13377- 815) + 273)e, 

+(5(13579 - 356977)77 - 21389)e2 + 107/(2047// - 8163) + 25434] , (315) 

(5?7(355r/ - 4843) + 7527)e? 



4. 3 i\E\e^i-et We,) 
am ■= 660(377- 1) 

+10r/(99377 + 6575) - 22522 , (316) 
25[52] -2zetWf^ + ^^i^-^ [(77(14577- 937) + 337)e? 

-22(r/(113r7 - 289) + 87)1 , (317) 
2~ 7|g|6^(77(62577 - 4849) + 1511)et Y^l 

"[^^1 ■= 99^^:33 ' (^^^^ 

z|i;|e^(577(12577 - 653) + 983)et>Vj, 
"[72] 44(3,7-1) ' ^^^^^ 



(4) 

= :TTA/rl^lV(l-3r/) <! >7a[fc2i-F[fe2i(M)!>, (320) 



6|g|e^(577(8r7-41) + 52) 

^ 55(3^^1^ ' (^'^^ 

6|i;|6^(577(877 - 41) + 52) 
«[iol ■= -6 + ^(^^ ' (322) 

^ , |£;|e^(57;(82r;- 1699) + 2711) 
"[20] .= -5 + 55(377 - 1) ' ^^2^) 

, 2 ^ . ^ |g|>VeV^(5^(2r/+163)-218) 
"[30] := -5et+5 + .... _^^ , (324) 



3|g|>V,V^(2077(777 - 128) + 831) 
"[40] fifi,.,oo ' (325) 
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^"[50] 


gli;! £2(207^(777 - 128) + 831) 

110(377-1) 


(326) 


(4) 


= 0, 


(327) 


(4) 
043 

o 


■srv, 16V^|i?l^/'We.<5™(l-2r;)r; 

_ „-i3KM „-i34>o V 

- " " 15e? 






I fc=2 fe=3 J 


(328) 






(329) 


'/3[33] 


:= 2 (e^ - 5e? + 4) , 


(330) 


^/3[43] 


:= i(-et + 37e?-36) , 


(331) 




:= i We' (336? -68) , 


(332) 


'/3[63] 


:= 70 We^^, 


(333) 


^/3[73] 


:= -30 Wf^, 


(334) 


^/^[33] 


:= -zWe,et(e?-4) , 


(335) 


^/3[43] 


:= 7We,et(l2-7e2) , 


(336) 


^^53] 


:= ^7W,,et(e^ + 3e?-4) , 


(337) 


/3[63] 




(338) 


^^[73] 


:= 30iet Wl^ , 


(339) 


(4) 
^42 


= 0. 


(340) 


(4) 
^41 


— e e X 

35et 






lfe=2 fc=3 J 


(341) 




•— ^ 7 






2 We' 

3 ' 


(343) 


A41] 


25 W^^ 
6 ' 


(344) 




2 ' 


(345) 


^/3[31] 


:= -iWe,et, 


(346) 


^/3[41] 


:= -^*We,et, 


(347) 


^/3[51] 


:= ^*e,We^,, 


(348) 


(4) 
^40 


= 0, 


(349) 
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(5) 
j55 


4 

r 9 9 >. 






\ Y ^"[fc5]-7^[fe5](tt) + Y ^^m^S[kb]{u) \ , 


(350) 




Ia:=0 fc=l J 




^"[05] 


:= I {et - 12e? + 16) , 


(351) 


^"[15] 


:= -I {et - \2e\ + 16) , 


(352) 


^"[25] 


:= -^j(73e4- 396e? + 400) , 


(353) 


3 

^[35] 


.— 6 I 


('354'! 


3 

a [45] 


:= --^J (et - 1) {et + 1) (83e^ - llSSe^ + 1840) , 


(355) 


3 

"[55] 


:= {Set - 1028e? + 2096) , 


(356) 


^[65] 


7 

• = — iyy6 (236? + 112) 


(357) 

you 1 ; 


^"[75] 


:= V<(3e?-8), 

U 


(358) 


^"[85] 


210iWi°, 


(359) 


^"[95] 


:= -70i Wif , 


(360) 




-5e^ + 20e^ - 16et 


(361) 


^a[15] 


We, 


^"[25] 


:= 0, 


(362) 


^Q;[35] 


et (-42e^ + 341e^ - 692e^ + 400) 

IZ, f Vet 


(363) 


^(5 [45] 


(et We J (I07ef - 564e2 + 560) 
12 


(364) 


^"[55] 


:= -let We,' (I9e4 - 1566? + 240) , 
o 


(365) 




^etWe,' (31e2- 176) , 


(366) 


^"[75] 


^etWe/(l5e? + 16) , 


(367) 


^"[85] 


:= MOetWe^^, 


(368) 


^"[95] 


-70etWeV, 


(369) 



(5) 

7^'' = 0, (370) 



-srM 8fy\Ef/'S^v{l-2rj)We, (e? - 4) 

j53 ~t3KM —23(t>a V "''^ ^ 

3e3 



\ 5Z ^"[fc3]-^[fe3] (w) + 51 ^'5[fe3]-7^S[fe3] (?i) > , (371) 
lfc=0 fc=l J 

"[03] (372) 
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^"[13] := (373) 

3 » (1736^404) 

"[^^1 60 (e?- 4) ' 

3 » (9ef - 25e^ + 16) 

"[^^1 ^= 30 (e^ - 4) ' 

. 7i (73ef - 301e2 + 228) , , 
- - 120 (eg- 4) ■ <™' 

120 (e?- 4) ' ^ ' 
3 77i>Vf 

"[63] g (g2 _ 4) ' (378) 

3 11* . ^ 

^5[23] := 0, (381) 

3- . (l26e^-495e^ + 404) 

"[^^1 ■= 60 (e?-4) ' (^^2) 

3~ . (59e^-183e^ + 124) 

20 (6^4) ' 

3 , e,Wf^ (896^ 284) 

"[^^1 ■= 40 (e?- 4) ' ^^^^^ 

3 22etW5 

(5) 

= 0, (387) 

fsl ^-.KM 16y#|^;|^/V<5,„(l - 2^) We. ' 

J — o C A 

3met 

fe^"[fci]-^[fci](") + II^"[fei]-^s[fci](M)| , (388) 

lfc=0 fe=l J 

^«[oi] := (389) 
3a[ii] := -i. (390) 

77j 

'«[2i] := (391) 
7iW2 

To" ' 

3 35iW2 

«[4il := (393) 
3a[5i] := (394) 



'«[3i] ^TT^. (392) 
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e* 

We, ' 


(395) 


^a[2i] 


0, 


(396) 


^Q;[31] 


7et (11 - 6e2) 

60 We, 


(397) 


Q![41] 


7 




3 ~ 
C*[51] 


7 , 

•- g^t'^'^etJ 




(5) 
j50 


= 0, 


(400) 


^55 


= 0, 


(401) 


(5) 


.-,4.^.„-.4..64v/il^l'^(5(^-i>' + i)., 

f ' ~ 1 






s N ^/3ri,/(i J^ri./ii fli) + N ^/3ri./ii J^cn./ii fit) / . 
lfc=2 fe=3 J 


(402) 


4« 
P[24] 


■ /' 4 o 2 1 


/ /I no \ 

(403) 


^/3[34] 


:= (-3e^ + 31e^ - 60e? + 32) , 


(404) 


''/3[44] 


:= -^»((e? + 3)(l5e?-34)e? + 76) , 


(405) 


'*/3[54] 


■■= ^*We' (3ef + 218e?-396) , 


(406) 


^/3[64] 


(I9e2 + 16) , 


(407) 


^/3[74] 


-^*We' (33e?-100) , 


(408) 


■*/3[84] 


:= -105iWi°, 


(409) 


P[y4J 


•= 35iyV^^ 


(410) 


P[34] 




(411) 


''/3[44] 


(-et WeJ(29e4-112e2 + 88) 
12 


(412) 


'^/3[54] 


:= -\ety^l (3e^-16e? + 8) , 


(413) 


■*/3[64] 


:= -fe,W.^^(e?-6), 


(414) 


''/3[74] 


:= -fe,Wl,(3e?-4) , 


(415) 


''/3[84] 


:= -70e*Wf^, 


(416) 


''/3[94] 




(417) 


(5) 


= 0, 


(418) 
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(5) 


,2;cA. .2.- 64yff|i?|V(5(r?-l)r/+l)W2 

45e2 

f ^ ^ ~ 1 






< 51 '*%2]-^[fc2] (m) + XI '^Afc2]-^S[fc2] (m) > , 


(419) 




I fc=2 fe=3 J 




'^/3[22] 


:= -i , 


(420) 


''/3[32] 


i^ (3e? - 10) , 


(421) 


4 o 


■■= -Y^*(3e?-8) , 


(422) 


^152] 


7 

• = — i fSp! - lie? + 8) 




P[62] 


3 ' 


(424) 


4« 
P[72] 




(425) 


''/3[32] 


et (e?-2) 


(426) 


'^/3[42] 


et (lle2 - 16) 


(427) 


"^52] 


et (3e^ - lie? + S) 


(428) 


''^62] 




(429) 




:= 5etW,5^, 


(430) 



(5) 



(5) 

S^^ = 0, (431) 



= yy^|i?|VWe,(5(r,-l)r? + l)etX |^^%o]-^5[feo]M| , (432) 

■■= -1, (433) 

^^42] := (434) 

%2] := (435) 



j66 _ g-j6/CA1 g-i60o V ' ^ 

3e? 



11 11 



''"[06] 


:= 4- 






''"[16] 


■■= -( 




1 


''"[26] 


60 




1 


''"[36] 


30 



'^"[fc6]-^[fc6] (■") + ''"[fc6]-^S[A:6] (u) \ , (436) 
. fc=0 fc=l J 

- 18e^ + 48e? - 32 , (437) 

e? - 2) {et - 16e2 + 16) , (438) 

(-437e^ + 3666ef - 7536e2 + 4384) , (439) 

(et - 1) (et + 1) (I5e^ - 444e^ + 1648e2 - 1408) , (440) 
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^"[56] 



''"[66] 



a[46] := (et-l)(et + l)(283e?- 7086e^ + 24208e?- 19680) , (441) 

Wl (I5e^ - 5598e^ + 334886^ - 37600) , (442) 

-^^l (719e^ - 19846^ - 12560) , (443) 

:= -^Wl (I605e^ + 8336e? - 44496) , (444) 

'"[86] := -^Wj°(369e? + 406) , (445) 

'"[96] := y Wif (51e?- 190) , (446) 

4«[i06] 924 (447) 

4a[n6] := -252 W^f, (448) 

a[i6] -2zet We, (36^ - 166? + 16) , (449) 



'a 

4 ;r 



5[26] 0, (450) 



4 ~ 



4 ;^ 



4 ~ 



4 ~ 



4 ~ 



-i(et We J (I20e? - 1251e^ + 3152e? - 2192) 



"[36] ^ , (451) 



ou 

i{et We J (329e? - 2823e^ + 5984e? - 3600) 



5[46] ^ ^ (452) 

4~ -^(et WfJ (55e^ - 501e^ + 1216e^ - 880) 
"[56] := ^ : (453) 

-z(et W^J (301e^ - 3296e? + 5360) 
"[66] := — ' 

z(et WJJ (27e^ + 704e^ - 1776) 
"[76] Y2 ' ^ ^ 

56zWe>t (2e^ + 51e^ - 53) 
"[86] := ^ — (456) 

4~ 28zW9et(l2e^-35e^ + 23) 
"[96] ^ , (457) 

4a[io6] := -672iWlfet, (458) 
4«[n6] 252ietWif, (459) 



(6) 

/^^ 0, (460) 



.-.4^A. ,-.4.n 64yi|r|£;|3^(5(^- 1)^+1) 



< X! ''^[fc4]-7^[/c4] (w) + X! '"[fc4]-^-S[M] (w) > , (461) 
lfe=0 fe=l J 

4a[04] := 8-8e? + e4, (462) 

^a[i4] := -6^ + 86? -8, (463) 

^a[24] ^ (-237et + 1096e? - 936) , (464) 

'"[34] := ^(-15e?-97ef + 456e?- 344) , (465) 
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"'a [44] 


• — 


"360 " ^""^ "'^ ^^^^^""^ ~ + 


12584) , 


(466) 


^"[54] 




^Wf^ (I65e4 _ 4264e? + 9384) , 




(467) 


4 

a [64] 




Wf, (31e2 - 356) , 

45 




(468) 


''"[74] 




iw« (57e?-170) , 
y 




(469) 


'^"[84] 


:= 


49W,\°, 




(470) 


''a [94] 


:= 


49 VV^^ 
3 ' 




(471) 


4 ~ 
"[14] 




-4zet Wet (e* - 2) , 




(472) 


^"[24] 




0, 




(473) 


^a[34] 




i{et We,) (40e4 - 217e2 + 234) 
15 




(474) 


^(5 [44] 


:= 


i{et We J (209e^ - 973e? + 874) 
45 




(475) 


'*Q![54] 




*(et WfJ (■^5e4 - 431e2 + 466) 




(476) 




oU 




4 ~ 




2i(et WfJ (251e2- 936) 
45 




(477) 


''a [74] 




KVe.et (3364 - 248j , 




/ A '-7Q\ 

(478) 


''a[84] 


:= 






(479) 


4 ~ 
Q![94] 




^^e W" 




y^ou ) 


(6) 
^63 


— 


0, 




(481) 


(6) 
j62 




■,rv, 32Jif|i?|V(5(77- 1)77+1) 








33e? 

r 7 7 >| 

lfe=0 fc=l J 


X 


(482) 


''"[02] 


I— 






(483) 


4 

«[12] 










4 

"[22] 




97 39e? 
30 20 ' 




(485) 


''"[32] 




^ {-blet + 95e? - 44) , 




(486) 


^"[42] 




-^(27e^-89e? + 62) , 




(487) 


^"[52] 




-^>Ve^.(21e?-26), 




(488) 


^"[62] 




175 We^^ 

36 ' 




(489) 



Copyright line will be provided by the publisher 



andp header will be provided by the publisher 



47 



4 

'^[72] 


25 Wf, 
12 ' 




('490'! 


4 - 
"[12] 


:= ~2iet We, , 




(491) 


''a[22] 


:= 0, 




(492) 


^"[32] 


let (406* — 137ej + 97) 

30 We, 




(493) 


4 ~ 
"[42] 


iet (I37e4 - 384e? + 247) 




(494) 


90 We, 




^"[52] 


z We,et (87e^ - 284e? + 


197) 


(495) 


DU 




4 ~ 
"[62] 


25. 




(496) 


4 ~ 
"[72] 


25 ,,,7 

- 12 




(497) 


(6) 
j61 


= 0, 




(498) 


(6) 


64 /lO^.^.g ^. 

= 33V^'^'^('^^-'^^ 


f ' 1 

+ 1) ^ "[fcO]-^[fcO](w) ) , 


(499) 


^"[00] 






(500) 


^"[10] 


-1, 




(501) 


^"[20] 


77 

~60 ' 




(502) 


'^"[30] 


7 We', 

10 ' 




(503) 


4 

"[40] 


35 W^ 
24 ' 




(504) 


''"[SO] 


7W4^ 

ft 




(505) 
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